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A FINITE ELEMENT FOR THERMAL STRESS ANALYSIS 
OF SHELLS OF REVOLUTION 


By Howard M. Adelman, Harold C. Lester, and James L. Rogers, Jr. 

Langley Research Center 

SUMMARY 

This report describes a new axisymmetric finite element for calculating static 
thermal stresses in general orthotropic thin shells of revolution. The element is geo- 
metrically exact and the thermal loading conditions allow for variations over the shell 
surface as well as through the shell thickness. 

The element is utilized in a number of sample calculations on a variety of shell 
shapes and loading conditions. Among the shells analyzed are two cylinders, a conical 
frustum, a truncated hemisphere, and an annular plate. Results from the present method 
are compared with either exact solutions which are developed in the report or with results 
from other methods based on finite differences. The results predicted by the present 
method were found to be in excellent agreement with the exact solutions and the finite- 
difference solutions. 


INTRODUCTION 

The capability to predict reliably the static stress of structural components sub- 
jected to thermal loads has become an urgent need for structural designers and analysts. 

Of particular interest are aerospace vehicles undergoing aerodynamic heating during 
atmospheric entry and aircraft structures heated by impinging jet engine exhaust. In many 
instances the aerospace vehicles of interest are thin shells of revolution. Closed-form 
analytical solutions to thermal stress problems for this class of structure have been 
obtained only for limited geometries and loading conditions (refs. 1 to 3). For more gen- 
eral geometries and loadings it is usually necessary to resort to approximate methods 
such as finite differences (refs. 4 to 6), numerical integration (ref. 7), or the finite- 
element method (ref. 8). Finite -difference methods and numerical integration techniques 
have been used with a good deal of success for thermal stress analysis of shells of revo- 
lution with a high degree of generality in the thermal loadings. 

The finite-element method has developed into the primary analysis tool for complex 
structures such as complete aerospace vehicles. Large comprehensive finite-element 
programs such as NASTRAN (ref. 8) have a variety of finite elements available to the user 
for modeling almost any conceivable configuration. However the elements themselves, 
especially the shell of revolution elements, are quite primitive both with regard to their 
shape and the generality of thermal loads that can be applied, particularly when compared 


with the generality available in the previously mentioned approaches. It appears then that 
the finite-element method can be made more useful with regard to thermal stress analysis 
by introducing new and/or improved elements which have the kind of generality that is 
available in the finite-difference and numerical integration approaches. Accordingly, an 
improved axisymmetric finite element for thermal stress analysis of shells of revolution 
has been developed. 

The purpose of the present report is to describe and evaluate the new finite element. 
The finite element is geometrically exact, a concept previously utilized for analysis of 
free vibrations (refs. 9 and 10) and of temperature distributions in shells of revolution 
(ref. 11). 

The report contains a description of the mathematical development of governing 
matrices for the element and describes a number of sample calculations carried out to 
verify the accuracy and versatility of the element. These calculations are performed for 
cylinders under two different loading conditions, a conical frustum, a truncated hemisphere, 
and an annular plate. Wherever possible analytical solutions (some derived herein) are 
used for checking results. There are three appendixes to this report. Appendix A con- 
tains a summary of the basic equations governing the finite element; appendix B contains 
the derivation of the thermal stress in an annular plate under a quadratically varying 
temperature; and appendix C contains a method for obtaining exact modal solutions for 
thermal stresses in freely supported cylinders. 


SYMBOLS 


a 


inner radius of annular plate (appendix B) 


a 0,k> a l,k> a 2,k’"l 
a 3,k> a 4,k> a 5,k J 


coefficients in polynomial displacement function for normal 
displacement w 


b 


outer radius of annular plate (appendix B) 


b 0,k> b l,k 1 

b 2,k> b 3,k J 


coefficients in polynomial displacement function for meridional 
displacement u 



(n = 0) 

(n/0) 


Cjj, Cj 2 > Cg 2 membrane stiffnesses 

Cgg in-plane shear stiffness 
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c 0,k> c l,k’^l 
c 2,k> c 3,k J 

D ll’ D 12’ D 22 

D 66 

E l> E 2 

e l’ e 2 

e 12 

J J t 
e l, 2’ e 12 

{F} 

< F k> 

(f } 

(f } 

G 

h 

[I] 

i = /3T 


coefficients in polynomial displacement function for circumfer- 
ential displacement v 

flexural stiffnesses 
torsional stiffness 

Young’s modulus for meridional and circumferential directions, 
respectively 

middle- surface strains in meridional and circumferential 
directions, respectively 

middle -surface shear strain 
total strains 

thermal force column matrix for complete shell 
element thermal force column matrix 

thermal force column matrix (see eqs. (12) and (13)) 

modal force column matrix (eq. (29)) 

shear modulus 

shell thickness 

identity matrix 


K number of finite elements used to represent shell 

K ll’ K 12’ K 22’ K 66 stiffnesses representing interaction between in-plane and out-of- 

plane strains 
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L 

M tl’ M t2 


M l» M 2’ M 12 


m 

N 


N ti, N t2 


n 

Rp R 2 


r 

[S ] 


s 

T 


T i’ 


T 

o 


meridional length of shell 

thermal moment resultants in meridional and circumferential 
directions, respectively 

moment resultants in meridional and circumferential directions 
and twisting moment, respectively 

meridional wave number for freely supported cylinder 

order of stiffness matrix and force vector after edge constraints 
have been applied 

thermal forces in meridional and circumferential directions, 
respectively 

circumferential wave number 
principal radii of curvature of shell 

radius of shell measured in plane normal to shell axis 
shell stiffness matrix 
meridional coordinate 
temperature 

temperatures at inner and outer fiber of shell, respectively 


[T^] matrix relating coefficients of assumed displacement shapes to 

degrees of freedom on edges of element 

T max maximum value of temperature 


T 1> T; 


stress resultants in meridional and circumferential directions, 
respectively 
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shear stress resultant 


T 12 

U 

u 

V 

v 

w 

[X] 


X 

{y> 


z 


a l» a 2 

P 

£ k 

e 

K V k 2 


K 
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[A] 


total potential energy 

meridional component of middle -surface displacement 
strain energy 

circumferential component of middle -surface displacement 

normal component of middle -surface displacement 

matrix which describes assumed form of variables appearing in 
strain energy 

meridional coordinate measured within single element 

column matrix containing unknown displacements, rotations, and 
derivatives thereof 

coordinate in direction normal to shell surface 

coefficients of linear thermal expansion in meridional and cir- 
cumferential directions, respectively 

rotation of shell generator, w' — — 

R 1 

length of kth finite element 
circumferential coordinate 

changes in curvature in meridional and circumferential direc- 
tions, respectively 

twist of middle surface 

diagonal matrix of eigenvalues of stiffness matrix 
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eigenvalue of stiffness matrix 


k 


k 


m 




{^ k > 


p 


a V a 2 


a 12 

w 

a 

Subscripts: 

k 

k+1 

Superscript: 

n 


•mode parameter for freely supported cylinder, hitt/L 
(appendix C) 

Poisson's ratio for meridional and circumferential directions, 
respectively 

column matrix whose elements are displacements and rotations 
and their derivatives at ends of kth element 

mass density 

meridional and circumferential stresses, respectively 

shear stress 

modal matrix for stiffness matrix 

potential energy associated with thermal loading 


number of element or number of first edge of element 
number of second edge of element 


circumferential harmonic number 


Primes denote differentiation with respect to s or x; T denotes transpose of 
a matrix. 


ANALYSIS METHOD 

In this section of the report, the governing matrix equations are derived. The 
procedure follows very closely the development in reference 9. The strain energy for a 
thin orthotropic shell of revolution is derived in terms of the three components of dis- 
placement of the shell middle surface and in terms of the shell temperature. A finite- 
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element representation of the shell is introduced in which the shell is represented by 
congruent slices of the shell and the displacement components are approximated by 
polynomials. The potential energy is minimized to obtain the governing matrix equations 
which involve two main contributions: the stiffness matrix and the thermal load vector. 
The stiffness matrix is identical to that presented in reference 10. The present develop- 
ment derives the thermal load vector. 

Derivation of Governing Equations 

Pertinent geometrical quantities are defined in figure 1. The components of dis- 
placement in the meridional, circumferential, and normal directions are denoted by u, v, 
and w, respectively. A point in the shell is defined by the meridional, circumferential, 
and normal coordinates s, 6 , and z, respectively. The two principal radii of curva- 
ture of the shell middle surface are and Rg and the distance from the axis of the 
shell to a point on the middle surface is r. 

The stress- strain relations for a linear orthotropic material including effects of 
temperature as given in reference 12 are 


E 1 t t E 1 T 
0-1 = (ej + M 2 e p " ( a 1 + 




1 - 


^2 


^2 


1 ^ ^2 2’ 


*2 = 


1 " ^2 


(e 2 + 


E 2 T 

( a 2 + ^i a i' 


> 


1 - Mj^ 2 


ct 12 " Ge 12 


(1) 


where the subscripts 1 and 2 refer to the meridional and circumferential directions, 
respectively, and = m 2 E^. 

The strain energy of the shell modified for a two-dimensional orthotropic material 
is conveniently written according to reference 13 as 



j J J (Vi 

Js J 6 J z L 


+ Cr, 


2 e 2 


T 12 e 12 


-( 


a j cr x 4* 


a 2 a 2 


)T 


r dz d6 ds 


( 2 ) 


It is assumed that lines originally normal to the middle surface remain straight, unex- 
tended, and normal after deformation. The total strains then have the following variation 
through the shell thickness: 
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( 3 ) 


e^(s,£,z) = e^(s, e) + zkj (s, 6) 
e 2 (s,#,z) = e 2 (s,0)+ z* 2 (s,£) 
e« 2 (s,£,z) = e 12 (s,0) + 2 zk 12 (s,0) 


Substituting equations (3) into equation (2), rearranging, and discarding a constant involv- 
ing the known temperature yield: 


u = v - n 


(4) 


where 


v = ||J ( C n e l + 2c l2 e l e 2 + C 22 e 2 + C 66 e 12 > r d0 ds 


+ ~ JJ ( D iri + 2D 12 K 1*2 + D 22*2 + D 66 ,< 12^ r d0 ds 

+ I JJ (K ll e l K l + K 12 (K l e 2 + K 2 e l } + K 22 e 2 K 2 + K 66 e 12*l2 )r dd ds (5) 

n = jj (N^ej + M tl Kj + N t2 e 2 + M t2 * 2 )r d& ds (6) 


In equation (5), V is the strain energy of the shell which is utilized in the development 
of the stiffness matrix in references 9 and 10 and presented in this report for complete- 
ness. The stiffnesses in equation (5) are needed in a later section of the report and are 
given in table 1. The rest of the present development consists of operating on the 
thermal potential energy 0 to obtain the thermal load vector. 

In equation (6), the quantities and are thermal forces and and 

M^ 2 are thermal moments which are defined as follows: 


E i ( a i + ^2 a 2^ 


% - f ^ 

k 1 ” 


T(s,0,z) dz 


(7a) 


M 2 
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(7b) 


c E„(a- + Mi a i) 

N t2 = I _l_i Li_ T(s, 0 .z) dz 

Jit 


MiM 2 


M t l 


M t2 = 


♦ E l( a l + ^2 a 2^ 
z !-^ 2 

E 2 (a 2 + /x l a l ) 


r 

Jz 1 ’ 


T(s,0,z) z dz 


T(s,0,z) z dz 


M 2 


(7c) 

(7d) 


The components of displacement and the temperature are assumed to have the following 
separable forms for each harmonic n: 


u(s,£) = u(s) cos nd 
v(s,0) ^ v(s) sin n 9 




> 


( 8 ) 


w(s,£) = w(s) cos n& 

T(s,#,z) = T(s,z) cos nd j 

The strain displacement relations used are those of Novozhilov as used in reference 9 and 
repeated here for completeness: 


Middle -surface strains: 


e. = u' ♦ I. 

R 1 


1 3v r ! 




w 


e 9 = -+ — u + . 

z r W r R 2 r 


e 1? = + v’ - 2- v 

u r M r 


(9) 
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Changes of curvature: 


K 1 = ~ 


w” + — u' 
Hi 


Ri 


r; 


k 2 = 


-_Li^ + J_^I-£lw’ + _£Lu 

r 2 30 2 rR 2 r rR l 


( 12 


-I Bw ' + r ’ ^. + _L^L + J_v' - _£l_ V 


r 30 p 2 30 rRj 30 R 2 


rRn 


(10) 


The thermal forces and moments as well as the strains and changes of curvature are 
written in separated forms consistent with equations (8) 

N tl = N tl cos nd 
N t 2 = cos n 0 

M u - M ti cos n0 
^2 = ^2 cos n0 

e l = e l cos 


> 


e 2 = e 2 cos n 0 


e 12 = e 12 s * n n6 


K 1 = *1 cos 


k 2 ~ k 2 cos 


k 12 = *12 sin 


( 11 ) 
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where N^, N^, Mq, and so forth are functions of s only. Substituting equations (9) 
and (10) along with equations (11) into equation (6) yields 

0 = C n j <^w, w', w”, u, u', v, v’) T {f} ds (12) 

where {f } is a 7 x 1 column matrix whose elements are given in the following equations: 


f l = iT N ti + ^ N t 2 + 7 M t2 
K 1 k 2 r 


f 2 = -r’M t2 


f 3 = -rM tl 


f 4 = r ’ N t2-^ R i M tl + M t2 

R^ 1 * 


f 5 = rN tl + — M tl 
R 1 


f 6 = rN t2 + — M t2 
R 2 


f 7 = 0 


(13) 


The integration with respect to 9 
made of the following results: 

C 2n 2 

C = cos nd de 

n Jn 


f 2n 2 

S n = I sim n& dd 
n L 


from 0 to 2 tt has been performed and use was 

C2tj (n = 0) 

1^77 (n > 0) 

r 0 (n = 0) 

t (n > 0) 
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The thermal forces and as well as the thermal moments and 

are dependent on the nature of the temperature variation through the shell thickness. (See 
eqs. (7).) One important case is that of a linear variation through the thickness. For 
later use, the thermal forces and moments will be evaluated for this special case. If the 
temperature varies linearly from at the inner face (z = -h/2) to T Q at the outer 
face (z = h/2), then 


T(s, z) = 


T o +T i 

2 



(14) 


where T q and T. are functions of s. Substituting equation (14) into equations (7) and 
integrating through the shell thickness give 


N tl - W C 11 + a 2 C 12^ 
N t2 = ( a l C 12 + a 2 C 22^ 


T o +T i 


T o +T i 


+ (Kn a i 

+ (k 12 ai 


+ K 12 a 2 ) 
+ K 22 a 2 ) 







(15) 


M tl " ^ a l D ll + a 2 D 12^ 


T - T. 


T + T. 


° 1 + ( K H a l + K 12 a 2> "° _i 


M t2 - ^ a l D 12 + a 2 D 22^ 


T_ - T s T_ + T, 

" 12^1 + ^ 22 ^ 


° 1 + tfC-a, + K 00 a 0 ) 12_1 


These forms for the thermal forces and moments complete the development of the column 
matrix {f } in equations (13) for the special case of a linear temperature distribution 
through the shell thickness. 

The next step in the development is to write Q for a single finite element in order 
to be able to arrive at the required thermal load vector for an element. The basis for 
these operations is described in references 9 and 10; for completeness, the essential con- 
cepts are summarized in appendix A. Combining equations (A2), (A5), and (12) yields 
the following form of the thermal potential energy for the kth element: 

°k = / k } T ( F k) <‘«> 
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where (F k > is the thermal load column matrix for the kth element as given by 


«V - c „ "k 1 


I 


: k/ 2 

£ k/ 2 


[X] T {f> dx 


(17) 


Formulation and Solution of Final Equations 

As in reference 10, the following conditions of compatibility are imposed at each 
element juncture: 


w k+l 

1 

w k + l 

Vi 


u k + l 

v k + l 


v k + l 

•< ^k+1 

► = « 

^k+1 > 

u k + l 


u k + l 

v k+l 


v k + l 

J^+l> 

k element 

A+k 


(k < K) 


(18) 


k+1 element 


Of the seven conditions of compatibility shown in equation (18), only the first four are 
strictly required to assure that the element converges. The last three conditions are 
imposed in order to improve the convergence characteristics of the element. Specifically, 
these conditions require continuity of strains and change of curvature at element junctures 
and evidence of the improved accuracy resulting from this continuity has been presented 
in reference 12. 

Two final observations are that the last three conditions are not correct for shells 
which have a step discontinuity in stiffness at an element juncture and that for problems 
in which the theoretical stress distribution has a slope discontinuity, the present element 
will display a "rounded corner" at the location of the slope discontinuity. 
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The total thermal potential energy is obtained by summing the element contributions; 

thus, 


a = 



n, 


(19) 


where is given in equation (16). When the summation indicated in equation (19) is 
carried out and use is made of equations (17) and (18), the thermal potential energy may 
be written as 

0 - (y > T iF> (20) 

where 

{F} thermal load column matrix of order 7(K + 1) 

{y} column matrix containing the unknown displacements, rotations, and deriva- 

tives thereof 

From reference 9, the strain energy may be written 

V = - (y} T [S] {y} (21) 

2 

where [S] represents the stiffness matrix which is a symmetric positive definite or 
positive semidefinite matrix of order 7(K + 1). The formation of the stiffness matrix has 
been described in references 9 and 10 and need not be discussed further. The formation 
of the thermal force column matrix for the complete shell from the element thermal force 
column matrices is accomplished by overlapping the last seven numbers in each column 
matrix with the first seven numbers in the adjacent (higher numbered) column matrix. 
Finally the column matrix {y} is the union of the column matrices as defined in 

equation (A6). 

The equations governing the behavior of the shell before edge constraints are 
applied are derived by minimizing the total potential energy. Thus, 


S(V-fi)=0 (22) 

The minimization in equation (22) is equivalent to the following matrix equation: 

[S ] {y} = {F } (23) 
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Edge constraints are incorporated by deleting appropriate rows and columns from [S] 
and rows from (F). In table 2, the rows and columns deleted for various edee conditions 
are given. It is noted that the degrees of freedom of the edges constitute the first seven 
rows of {y} and the last seven rows of {y} . Further, it is only the first four of these 
in each group namely w, u, v, and /3 that are involved in the definition of edge con- 
straints. Accordingly, as shown in table 2 the only rows and columns that are ever deleted 
from the [S] matrix and the (F> matrix are among the first four and 7K +1, 7K + 2, 
7K + 3, and 7K + 4 rows and columns. 

The stiffness matrix is nonsingular if and only if the edge constraints are sufficient 
to prevent all possible rigid body motions of the shell. In this case, equation (23) may be 
solved by any of a number of standard methods for solving sets of simultaneous linear 
equations. Unfortunately in practical shell problems there are situations wherein the edge 
conditions are not sufficient to prevent all rigid body motion and the standard methods are 
not applicable. For example a free-free shell under a self-equilibrating set of loads has a 
singular stiffness matrix and more importantly (for this paper) so does a shell with an 
insufficient set of edge constraints under any applied thermal loading. Equation (23) cannot 
be solved for either of these classes of problems without a solution technique which is not 
dependent on the nonsingularity of the stiffness matrix. For this reason the solution of 
equation (23) is carried out by a method normally used for dynamics problems — that is, 
an eigenvector expansion. Let 

(y) = [0] {q} (24) 

where W is the matrix of eigenvectors of the stiffness matrix in which each column is 
an eigenvector denoted by {<£} and {q> is the column matrix of modal coordinates. 

Thus, {<£) and A satisfy the equation 

[SK<£}=\{0} (25) 

The eigenvectors are normalized so that 

[<£] T W = [I] (26) 


Then 


[<£] T [S] [0] = [A] 


(27) 
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where 


X 


1 


^2 


[A] = 


(28) 


k 


N 


T 

Substituting equation (24) into equation (23), pre multiplying by [<£] , and making use of 
equation (27) give 


[A] (q> = (f } (29) 

_ 'T' 

where (f } = W {F> is the generalized force column matrix in which the ith element of 
(f } is the component of force that would deform the shell into the shape characterized by 
the ith eigenvector of [ S ] . The solution of equation (29) is as follows 

\h!\ <v o n 

<li = \ f (i = 1, 2, . . N) (30) 

l0 (Aj=0)J 

The column matrix (y> is then obtained from equation (24). 


Displacement and Stress Recovery 


The coefficients of the polynomials representing the displacement field are obtained 
from equation (A5). The displacement components u, v, and w are then obtained by 
use of the equations: 


w = 
u = 

V = 


a 0,k +a l,k x + a 2,k x2 + a 3,k x3 

b 0,k + b l,k x + b 2,k x2 + b 3,k x3 

c 0 > k +c l f ^ +c 2 > k^ + c 3,k‘ 3 


4 5 

+ a 4,k X +a 5,k X 

> 


J 


(31) 
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The middle -surface strains and changes of curvature are obtained from equations (9) and 
(10). The fiber strains are obtained from equations (3) by substituting ±h/2 for z and 
are as follows: 


4(m>±|)= e^s,#) '< 1 (s,0) 

*2 ^s,<9,± — ^ ”■ ©2 (®> $) ~ K2(s,e) 

3 12( s,0 » ± ^ = e 12^ s ’^ * ^ 




(32) 


The plus (+) sign denotes the outer fiber and the minus (-) sign denotes the inner fiber. 
The stress resultants Tj, T 2 , and Tj 2 and the moments resultants Mj, M 2 , and 
Mjg are defined as follows: 


T l = 


T 2 = 


T , 0 = 


fh/2 

a l 

dz 

m i= | 

r h/2 a.z dz 

'-h/2 


j 

-h/2 

f h /2 

2 

dz 

m 2 = 

r h/ 2 

o-gZ dz 

'-h/2 



-h/2 

rh/2 

'-h/2 

12 

M 12 = 

fh/2 

^l 22 

■'-h/2 


(33) 


Substituting equations (1) along with equations (3) into equations (33) and using the definition 
in table 1 yields the following expressions: 
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T 1 ~ C ll e l + C 12 e 2 + K H k 1 + K 12*2 “ N tl 


T 2 " C 12 e l + C 22 e 2 + K 12*l + K 22*2 " N t2 


T 12 C 66 e 12 + K 66*12 


l l " 11*1 + 12*2 “ll e l + 12 2 “ "tl 


M 2 ~ D 12 k 1 + D 22*2 + K 12 e l + K 22 e 2 " M t2 


M 12 ~ \ D 66*12 + \ D 66 e 12 


The displacements and stress and moment resultants are evaluated for specified values 
9 . according to the following formulas: 


u(^,s) = u(s) cos n£j 


v(<9.,s) = L v(s) sin n0. 


w(<9.,s) = w(s) cos ne^ 


(i = 1, 2, . . .) (35c) 


T i v> = E 


T.,(s) cos ne. 

1 l 


(i = 1, 2, . . .) (35d) 


T 2 (6>i,s) = L TgCs) cos n^ 


(i = 1, 2, . . .) (35e) 


T 12^i’ s ^ " L T^s) sin n6L 
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<» - 1 , 2 , . . .) 


M i< v> = E 

n 


M^(s) cos n#,. 


(35g) 


M 2 (£.,s) = 2 ^ M 2 (s) cos n£. (i = 1, 2 , . . .) (35h) 

n 

m 12 ( V> - 22 M 12 (s) sin (i = 1, 2, . . .) (35i) 

n 

The range of summation on n is over the number of circumferential harmonics 
considered. 

EVALUATION OF THE ELEMENT 

In order to verify the validity and versatility of the new element, sample calculations 
were performed for a variety of shell geometries and thermal loads. Analytical solutions 
were available for some of the models. In the following sections are presented descrip- 
tions of the shells, the methods of evaluating numerical results, the finite -element repre- 
sentations, and detailed evaluations of the results. 


Description of Shells Analyzed 

The geometry of the shells analyzed are illustrated in figure 2 along with the cor- 
responding temperature loading functions and the shell material properties. For cylin- 
ders I and II and the truncated hemisphere, the temperature is constant through the shell 
thickness. For the conical frustum and the annular plate, the temperature varies through 
the thickness. For all shells except cylinder II, which has freely supported edges, the 
edge conditions are clamped-clamped. Also, the temperatures are axisymmetric except 
for cylinder II. A brief description of each shell is given as follows: 

Cylinder I is a cylindrical shell having a temperature distribution with a discontinu- 
ous slope at the midlength of the cylinder (s/L = 1/2). This discontinuity occurs when the 
midlength circumference is maintained at a prescribed temperature. The temperature dis- 
tribution (derived in refs. 3 and 11) illustrates the behavior of the present finite element 
for a stress distribution having a slope discontinuity. 

Cylinder II is a cylindrical shell with a temperature distribution which is constant 
along the length of the shell but varies around the circumference. This temperature dis- 
tribution (derived in refs. 3 and 11) results when a generator of the cylinder is maintained 
at a constant temperature. This example was chosen to illustrate the application of the 
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finite element to a nonaxisymmetric temperature distribution and also to illustrate an 
application where, because of the freely supported edge conditions, the stiffness matrix is 
singular. 

An orthotropic 30° conical frustum has a temperature distribution which varies 
quadratically along the shell length and linearly through the thickness. This distribution 
is an approximation to that in the vicinity of the nose of a missile undergoing aerodynamic 
heating (ref. 14). This example illustrates the use of the present element for orthotropic 
shells, for a nonconstant shell radius, and for an applied thermal bending moment. 

A truncated hemisphere having a temperature distribution (derived in ref. 11) cor- 
responding to a uniform heat flux over the outer surface of the shell was chosen to illus- 
trate a calculation for a shell with a curved generator such that both principal radii of 
curvature are nonzero. 

An annular plate with a temperature varying linearly through the thickness and 
quadratically from the inner radius to the outer radius illustrates the application of the 
element when thermal bending moments are applied as the loads and for the limiting case 
of an annular plate. 


Method of Evaluation 

The evaluation of the element was accomplished by using a computer program to 
calculate stress and moment resultants and then by comparing the results with other 
solutions for the same configurations. Check solutions were obtained from the following 
sources: 

For cylinder I, an analytical solution was available from reference 3 

For cylinder II, a converged modal solution was generated based on the procedure 
described in appendix C, and results were obtained for harmonics n = 0 to 4 

For the orthotropic conical frustum, check results were available from the finite- 
difference computer programs of references 5 and 6 

For the truncated hemisphere, a solution was available from reference 5 

For the annular plate, an analytical solution derived in appendix B was utilized 

Detailed comparisons between results with the present method and results from these 
sources are presented in a later section. 

Finite- Element Representation 

The number and spacing of elements in the shell representations are given in table 3. 
These representations were arrived at by the following procedure: first, each shell 
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was analyzed with a representation of 10 equally spaced elements and results were com- 
pared with the appropriate check solution. If satisfactory agreement was obtained (as 
for cylinder II and the annular plate), no refinement was performed. If the results were 
not satisfactory, a more refined representation of the shell was made in the vicinity of the 
edges of the shells where the steepest gradients in the stress distributions are found. The 
refinements were similar to those done in the modal stress study of reference 12. The 
final representations used to give the results in the present report are those in table 3. 

Evaluation and Discussion 

Calculated results from the present finite -element and other analyses are presented 
in tables 4 to 8 and in figures 3 to 8. In this section, an evaluation of the performance of 
the element for each set of calculations is presented. Before proceeding to discussions 
of individual shells, two general observations are made. 

(1) For all shells except cylinder II, results are axisymmetric; that is, they corre- 
spond to n = 0 because of the axisymmetric character of the temperature loads 

(2) For all shells except the annular plate, the moment resultants are negligibly 
small and are not presented 

Cylinder I .- Nondimensional meridional and circumferential stress resultants pre- 
dicted by the present method for cylinder I are presented along with analytical results 
from reference 3 in table 4. The stress resultant distributions are necessarily symmetric 
about the midlength of the cylinder as a consequence of the symmetry of the load and the 
edge conditions. The analytical solution for this example as given in reference 3 shows a 
constant meridional stress resultant T^, whereas the present finite -element solution 
exhibits a slight oscillation about the constant value. The maximum error of 0.6 percent 
occurs at the midlength of the cylinder. The comparison for the circumferential stress 
resultants T 2 also shows close agreement, the largest error being 0.4 percent. In 
figure 3, the stress resultants are plotted with the analytical solution being represented by 
the solid curves and the present finite-element solution indicated by the circles. It is of 
particular interest to observe that the analytical solution for T 2 has a cusp — that is, a 
discontinuity in the slope at the cylinder midlength as a result of the discontinuity in the 
slope of the temperature distribution at that location. Since the displacement field assumed 
for the finite element implies continuity in the slope of the stress resultants it might be 
expected that the element would have difficulty converging in the vicinity of such a cusp. 
However as indicated in figure 3, no such difficulty was found, and although there is a local 
rounding of the cusp, this rounding was too negligible to be plotted and had no appreciable 
effect on the accuracy of the stress resultant. 

Cylinder II .- Circumferential stress resultants T 2 for cylinder n for values of n 
from 0 to 4 are shown in table 5. In this table, the exact solution vas obtained by the modal 
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superposition method described in appendix C. The stress resultants from the present 
method are quite close Jo the exact solution for all values of n in the calculations. 

Further there is only a slight loss of accuracy with increasing values of n. The results 
are also shown in figure 4 wherein the present results are indicated by the symbols and 
the exact solution shown by the solid curves. Because the resultants are symmetric about 
the midlength of the cylinder, results are plotted for only half the shell. It may be observed 
from figure 4 that the harmonics of circumferential stress resultants decrease in magni- 
tude toward the midlength of the cylinder and the rapidity of the decay appears to increase 
with n. A final observation from figure 4 is that all of the harmonics pass through zero 
at about s/L = 0.125 for this particular cylinder. 

Figure 5 illustrates the results of performing the Fourier summation of the har- 
monics by using equations (35) at values of 6 of 0°, 45°, and 90° where the angle is meas- 
ured from the heated generator. Results in this figure are from the present method only. 

The meridional stress resultant Tj was found to be less than 1 percent of the 
circumferential resultant for all harmonics. This occurrence may be explained by the 
fact that the freely supported edge conditions require T^ to vanish at each edge, and for 
short cylinders such as that used in the present example, Tj never develops a significant 
value. For larger cylinders with freely supported edges, Tj is not small and this point 
is discussed further in appendix C. 

Orthotropic 30° conical frustum .- The meridional stress resultant T^ for the 
orthotropic 30° conical frustum which was computed by the present method and solutions 
from the methods of references 5 and 6 are presented in table 6. The agreement between 
the present method and reference 6 is excellent with a maximum difference of less than 
1 percent which occurs at the smaller edge of the shell. The difference between the 
present method and reference 5 is somewhat larger, being about 5 percent at the smaller 
edge and about 3 percent at the larger edge. The solutions are plotted in figure 6 where 
it can be seen that all three results are in close agreement except in the immediate vicinity 
of the edges of the shell. 

The circumferential stress resultant was negligibly small away from the edges of 
the shell and significant values occurred only within 1 percent of the shell length from 
the small edge and 3 percent of the shell length from the large edge. These edge results 
are summarized in the following table: 


Location 

Dimensionless circumferential stress resultant 
in a conical shell, T 2 /E ha T max 

Present 

analysis 

SADAOS 
(ref. 6) 

SALORS 
(ref. 5) 

Smaller end 

-28.17 

-29.44 


-28.89 

Larger end 

-12.24 

-10.82 


-12.27 
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The results are in fair agreement, with the present solution lower than either reference 5 
or 6 at the smaller end but bounded by those results at the larger end. The maximum dif- 
ference is about 13 percent with respect to the results of reference 6 at the larger end. 

One interesting and unexpected result in this calculation is that even though there is a 
temperature gradient through the shell thickness and, consequently, a thermal moment 
applied, the moment resultants were substantially zero. This result was obtained by all 
three analyses and cannot be rationalized by the authors. 

Truncated hemisphere .- Meridional and circumferential stress resultants for the 
truncated hemisphere which were computed by the present formulation and those from the 
method of reference 5 are presented in table 7. The agreement is generally good. The 
only significant differences between the two results occur in the values of the meridional 
stress resultant at the edges of the shell. There is a difference of about 5 percent at the 
larger diameter (s/L = 0) and about 10 percent at the smaller diameter (s/L = 1). There 
is good agreement between the two solutions for the circumferential stress resultant with 
a maximum difference of 2 percent which occurs at the larger edge. 

f 

The results are shown graphically in figure 7 where the prfesent solution is repre- 
sented by the symbols and the results of reference 5 are shown as the solid curves. An 
interesting and encouraging aspect of the comparison of results for this case is that, for 
the dominant stress resultant T 2 , there was excellent agreement and the significant 
differences occurred in the smaller stress resultant T^ at a location where it was 
substantially zero. 

Annular plate .- Meridional and circumferential moment results for the annular 
plate from the finite-element calculation are presented in table 8 together with the 
analytical solution derived in appendix B. For the meridional moment resultant compari- 
son, the solutions agree to within 0.3 percent except at s/L = 0.2 where the present 
method is 1.2 percent lower than the analytical solution. For the circumferential moment 
resultant comparison, the solutions agree to within 1 percent except at s/L = 0.2 where 
the two-place agreement at the inner edge and four-place agreement at the present solution 
is 4 percent lower than the analytical solution. 

CONCLUDING REMARKS 

A new finite element for thermal stress analysis of orthotropic thin shells of revolu- 
tion is described. The element is geometrically exact in that no approximation of the 
shell generator is required. Further, provision is made for the applied thermal loads to 
vary over the shell surface and through the shell thickness in any mathematically 
describable fashion. 
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Sample calculations are carried out on a variety of shell models and results from 
the present method are compared either with analytical solutions or with numerical results 
from independent analyses. Among the shell configurations analyzed are cylinders under 
two different thermal loading conditions, a conical frustum, a truncated hemisphere, and 
an annular plate. For these calculations the results from the new finite element generally 
are in good agreement with the check solutions. The only significant inaccuracies were 
observed in some stress resultants near the clamped edges of the conical frustum and the 
hemisphere where the stress gradients were particularly severe and consequently where 
finite-element analyses usually require significant refinement before convergence is 
obtained. 

The nature of the finite element is such that continuity of strains and curvatures is 
maintained at all points along the shell. Consequently in a shell where there is theoreti- 
cally a slope discontinuity or cusp in a stress or moment resultant the results predicted 
by the present method will show a slight local smoothing of the cusp but this has no 
adverse effect on the accuracy of the results. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., June 30, 1973. 
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APPENDIX A 


REPRESENTATION OF A SHELL OF REVOLUTION BY GEOMETRICALLY 

EXACT FINITE ELEMENTS 

The basic assumptions underlying the finite element and a detailed derivation of the 
stiffness matrix are given in references 9 and 10. The purpose of this appendix is to 
summarize the concepts, conventions, and equations needed in the development of the 
thermal load vector for a geometrically exact finite element. 

The following notation is used (fig. 9): 

K total number of elements used to represent shell 

e k meridional length of kth finite element 

x coordinate in kth element, measured from center of element so that 



A subscript notation is established in which a subscript k, when used with a quantity 
such as u, v, w, p , R^ or their derivatives, implies that the quantity is evaluated 

at the first edge of the element x = - -g-. Similarly, a subscript k+1, when used with 
such quantities, implies that they are evaluated at the second edge of the element 

x = ~ 2 ' In all other uses of the subscript k, for example, for matrices and vectors, it 

means that the quantity with the subscript is evaluated for the kth element. The displace- 
ment components or derivatives thereof within the kth element are approximated by poly- 
nomials as follows (ref. 10): 




w’ 


w" 
-< U 
u' 


> = [x] { V 


V 

V’ 

V y 


(A2) 
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where 


[xj = o 


{y} T = 


X 

*2 

x3 

x4 

X 5 

0 

0 

0 

0 

0 

0 

0 

1 

2x 

3x2 

4x 3 

5x 4 

0 

0 

0 

0 

0 

0 

0 

0 

2 

6x 

12x 2 

20x 3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

X 

X 2 

X 3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

2x 

3x 2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

X 

X 2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

2x 

,k’ a 

l,k ,a 

2.k ,a 

3,k’ a 4,k’ a 5,k’ 

b 0,k 

,b l,k’ 

b 2,k’ 

b 3,k’ 

c 0,k’ 

c l,k’ 

c 2,k’ 


The coefficients of the polynomials are related to the variables at the edges of the element 
by the following relation: 

{ V = [T k ]{ ^k } (A5) 


where the matrix [T^l is given in table 9 and 


(f k > = 


L w k’V T lt’ fi k’ u k’ T k>^k’ w k+l’ u k+:’ v k+l> /, k+l’ u k+l’ v k+l-' 3 lt+X_ 


The reader will notice that two new quantities have been introduced into (<f }. These are 
the meridional rotation fi and its derivative p' (which is identical to the change in 
curvature The reason for the introduction of these quantities in place of w' and 

w" is that in order to perform the necessary superposition of element matrices and vec- 
tors, a set of compatibility conditions at element junctures must be formulated and these 
conditions are convenient to express and are physically meaningful in terms of p and 
P' rather than w* and w". The matter of the form of the compatibility conditions at 
the element junctures is discussed further in the main body of this paper. 
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APPENDIX B 


ANALYTICAL SOLUTION FOR THERMAL STRESSES IN AN 
ANNULAR PLATE UNDER PURE BENDING 


The governing differential equation for the plate is derived in accordance with the 
following assumptions: 

(1) The plate is isotropic, thus 


a l " 




E 


1 " 



^i = 


= M 




(Bl) 


(2) The temperature loading is axisymmetric 

(3) The in- plane displacement components u and v are neglected 

For an annular plate having an inner radius a and an outer radius b, the functions 
describing the reference surface are 


r = a + s 


r' 


= 1 




R 1 



(B2) 


The principle of minimum potential energy is used to obtain the governing equation. The 
statement of the principle is 


S (V - n ) = 0 


(B3) 


Based on assumptions (1), (2), and (3) and equations (B2), equation (B3) may be written 


8 



D 



2j± 

r 


w'w” + 



+ 


2M t 

D 



dr = 0 


(B4) 
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where 


D - Eh3 
12(1 - M 2 ) 

and the prime indicates a derivative with respect to r. Carrying out the variations as 
indicated in equation (B4) yields the required differential equation 


w 


IV 




- rM^' 


Dr 


(B5) 


Letting the right side of equation (B5) be represented by a function M(r) and regrouping 
the left side yield 


I A 

r dr 


A I A 

dr r dr 



M(r) 


(B6) 


A solution to equation (B6) is sought which corresponds to a constant value of M. The 
general solution to equation (B6) with M equal to a constant is 


w = A + Br^ + cr^ log r + d log r + 


64 


(B7) 


The changes in curvature and the moment resultants are then calculated by use of equa- 
tions (10) and (37), respectively. The boundary conditions are taken to be those of 
clamped edges 


w(a) = w’(a) = 0 
w(b) = w'(b) = 0 


(B8) 


Substituting equation (B7) into equations (B8) gives the following equation for determining 
the constants: 
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1 

2 

a 

a 2 log a 

log a 


^ *N 

A 


r 

Tf< 

s ® 
1 

V 

0 

2a 

a(l + 2 log a) 

1/a 


B 

► 

_ Ma 3 
16 

> 

1 

b 2 

b 2 log b 

log b 

l 

c 

1 

4 

_ Mb^ 
64 

0 

2b 

b(l + 2 log b) 

1/b 


d 


_ Mb 3 

16 J 




- 





(B9) 


Return now to the determination of a temperature distribution corresponding to a constant 
value of M and set 


-Ml - rMl' 

L = M = Constant 

Dr 


then the differential equation for may be written as 


iL(rM{) = -MDr 

Integrating and discarding the constants of integration give 


(BIO) 


M t = 



(Bll) 


Using equation (Bll) along with the expression for the thermal moment in terms of tem- 
perature (eqs. (15)) gives 


M t = ^(1 + M HT 0 -T.)=- 


MD o 
r £ 

4 


(B12) 
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Thus a temperature distribution corresponding to a constant value of M is one in which 
the temperature difference between the inner and outer fiber of the plate varies quadrati- 
cally with r. Letting 



results in 


and 


M = -4a(l +m)£ 
h 


M _ aD(l + M )/^r 2 

t 1 


(B13) 


(B14) 


(B15) 
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EXACT MODAL SOLUTIONS FOR THERMAL STRESSES IN 
FREELY SUPPORTED CYLINDERS 

In this appendix the nonaxisymmetric state of thermal stress in a freely supported 
cylindrical shell is determined by a superposition of the free vibration modes. The 
potential energy for the shell is expressed in terms of middle -surface displacements and 
the temperature distribution in the shell. The temperature is assumed to be constant 
through the thickness. The displacement components and temperature distribution are 
represented in terms of their circumferential Fourier harmonics and each harmonic of 
stress is determined independently of the others in terms of the same harmonic of 
temperature. 


Symbols 

Symbols not previously defined in the main text and/or those unique to this appendix 
are defined as follows: 

summation index, k = 1, 2, 3 (see eqs. (C4)) 
stiffness matrix, 3x3 

stiffness coefficients (see eqs. (C7)) 

number of meridional half-waves 
mass matrix 
generalized force 

generalized coordinates (see eqs. (C4)) 

_ TO (n = 0) 

Sn ~ U (n/0) 

t time 

T* 

kinetic energy 


k 

[R] 


K ll> K 12’ K 13’ 
*22' *23, *33 


m 

CM] 

P 

nm 

%mk 
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a nmk’ ^nmk’ y nmk 

{a} 

"nmk 

Subscript: 

n 


modal amplitudes (see eqs. (C4)) 
column matrix of modal amplitudes 
natural frequency of nmk vibration mode 


nth circumferential harmonic 


General Method 


The strain-displacement relations used in this appendix are those of Kraus 
(ref. 15, p. 32) and differ from those in the main text (eqs. (10)) only in the k ^ term - 
The expressions used are 



zw" 




e 


t 

2 


I *1 + 2 + z/J- ll 

r 26 r l r 2 'de 





2H. + v' 

dd 


2 3w’ _ 1 V A 
r r / 


j 


(Cl) 


The strain energy U is derived in a manner analogous to that of the main text and is 
given by 


U =Z! U n 

n=0 


(C2) 


where U n is the contribution to the strain energy due to the nth circumferential harmonic 
of displacement and temperature and is given by the following equation: 
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where 

f 0 (n = 0) 

S nH 

|^7T (n ± 0) 

In equation (C3) the first integral represents the membrane strain energy, the second 
integral represents the flexural strain energy, and the third integral represents the 
thermal strain energy. 

The displacement components u n (s), v n ( s )> anc * w n^ s ^ are ex P an< ^ e< ^ i n terms of 
the free vibration modes of a cylinder with freely supported edges (ref. 9). This expansion 
is as follows: 


M 3 

u n (s) = 11 L a nmk cos ^ W 


m=l k=l 


(C4a) 
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M 3 


v n (s)= Z Z 


n • m77S „ 

^nmk 3 ; %mk 


(C4b) 


m=l k=l 


M 3 

W n (s) = Z Z 7 nmk sin q nmk 

m=l k=l h 


(C4c) 


Here m represents the number of meridional half-waves in the mode shape and M is 
the largest value of m in the expansion. The inner summation with index k is used to 
account for the fact that corresponding to a given nodal pattern defined by n and m 
there are three modes. The relative amplitudes of u n , v n , and w n in a given mode 
are given by a nmk , /? nmk , and ? nmk , respectively, and the generalized coordinate 
associated with a mode is denoted by <l nmk * The coordinates <l nrnk are calculated by 
the principle of minimum potential energy. Substituting equations (C4) into the strain 
energy in equation (C3), integrating with respect to s, and minimizing the resulting 
expression with respect to each of the unknowns q nmk give the following formula: 


^nmk 


Eahr 
1 - * 


P |-\ a 
nm 1 ~ ’ 


'nmk \ 


m a nmk ~ nmk ~ J 
r r / 


K ll a nmk +K 22 /3 nmk +K 33 7 nmk +2K 12 a nmk /3 nmk +2K 13 a nmk :K nmk + 2K 23 /3 nmk 7 nmk 


where 


(C5) 


m 


mrr 


(C6) 


nm 


r L 

= T n (s) sin \ m s ds 


and the stiffness coefficients are defined as follows: 


C Lr / o , 2 

K 11 = (Ck„ + c Lul — 
11 2 \ m 2 2 


(C 7a) 
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K 33 = 


C n Lr 


— +D^ + BIL + 2 m dIL- X^ + 2D(1 - m)— X^ 
2 m 4 2 m 2 m 

r r r r 


(C7f) 


where P mn is the generalized force for the mode characterized by n and m. The 
modal amplitude.s a nmk , /3 nmk , and 7 nmk are all that now remain to be calculated. 
They are computed by a modal analysis described in the next section. 

Determination of a nmk , /3 nmk , and / nmk 

The relative amplitudes a nmk , /J nm k’ and r nmk are obtained by a normal mode 
analysis based on the shapes shown in equation (C4). Making use of these equations, but 
noting that for free vibrations the modal coordinates <l nmk instead of being constant 
have simple-harmonic time variations, gives: 


M 3 

U (s,t) = / / a COS X s e 

n v ’ ' nmk m 

m=l k=l 
M 3 

v n (s ' l) ' £ Z ^nmk sin 6 
m=l k=l 

M 3 

w n (s ’ t > = £ II y nrak sin Xn s e 
m=I k=l 


i ^ 1 1 
nmk 


‘“nmk 1 > 


^ nmk 1 


(C8) 
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where ^ nmk is the natural frequency associated with the 
the triplet n,m,k. The strain energy U and the kinetic 


mode shape characterized by 
energy T* are written as 



T — 2 - V| tyi Irt 

{£} [K]{a}e nmK 


(C9) 


T* 
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2 

“nmk 


{a} T 


_ 2 &> 

[M]{a}e 


nmk*- 


(CIO) 


where 


{a} = 1 


nmk 


nmk 


nmk J 


and 

[M] = fI] 

2 


with [I] being a 3 x 3 identity matrix. Applying the principle of minimum potential 
energy along with equations (C9) and (CIO) gives the following modal equation: 


[K] {a} 


2 

"nmk 


[M] {a} 


This eigenvalue problem is solved by use of the threshold Jacobi method described in 
reference 9. 


Effect of Length on Thermal Stresses in a Cylinder 
With a Heated Generator 

In reference 3 a solution is obtained for the axial stress resultant in a cylinder of 
infinite length and used for predicting such a stress in finite cylinders. It was of interest 
in connection with the present work to assess the range of cylinder geometry for which the 
infinite cylinder solution is applicable to finite cylinders. Thus, a number of calculations 
using the exact modal method of appendix B were carried out in which all the shell 
dimensions were fixed except the length -to -radius ratio which was varied from 1 to 1000. 
The results are shown in figure 10. 
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The infinite cylinder solution is denoted by the broken line. It is of course con- 
stant. For L/r = 1, the axial stress is 0. For L/r = 10, the curve attains a compres- 
sive stress at the center of the cylinder which exceeds the infinite cylinder stress by 
about 9 percent. For L/r = 100 and 1000, the stress is constant and equal to the infinite 
cylinder solution outside a narrow boundary layer adjacent to each edge. At the beginning 
of boundary layer, a value of stress is reached which exceeds the infinite -cylinder solution 
by about 4 percent. The conclusion from this study is that the infinite -cylinder solution is 
a good approximation for the stresses in the interior of a long cylinder (L/r greater than 
or equal to about 100) but underestimates the stress near the edges. The infinite cylinder 
solution is furthermore found to be completely erroneous for cylinders having L/r less 
than 10. 
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TABLE 1.- STIFFNESSES OF AN ORTHOTROPIC SHELL 
AS USED IN EQUATION (5) 

[Subscripts 1 and 2 refer to meridional and 
circumferential directions, respectively] 



E 


1 



M 1E 2 

1 - Mi Mo 


1 - M-i Mo 


dz 


dz 




TABLE 2.- EDGE CONSTRAINTS 
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TABLE 3.- FINITE -ELEMENT REPRESENTATIONS OF SHELLS 
USED FOR SAMPLE CALCULATIONS 



Nondime ns ional element length. 

e k /L, for - 


Cylinder I 

Cylinder II 

Conical frustum 

Hemisphere 

Annular plate 

14 elements 

10 elements 

19 elements 

14 elements 

10 elements 

0.025 

0.1 


■ 

0.1 

.025 

.1 


■ 

.1 

.025 

.1 



.1 

.025 

.1 



.1 

.050 

.1 

.0166 

.12 

.1 

.050 

.1 

.0166 

.12 

.1 

.300 

.1 

.0166 

.12 

.1 

.300 

.1 

.1666 

.12 

.1 

.050 

.1 

.1666 

.12 

.1 

.050 

.1 

.1666 

.12 

.1 

.025 


.1666 



.025 


.1666 



.025 


.0166 



.025 


.0166 

■ 




.0166 





■ 





■ 





1 
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TABLE 4.- STRESS RESULTANTS FOR A CLAMPED-CLAMPED 
CYLINDER WITH A HEATED MIDLENGTH 


s/L 

T i / 

EhaT 

Tg^Ehat 

Present 

solution 

Analytical 
solution 
(ref. 3) 

Present 

solution 

Analytical 
solution 
(ref. 3) 

0 

-0.548395 

-0.548412 

-0.137099 

-0.137103 

.1 

-.548250 

-.548412 

-.235410 

-.235451 

.2 

-.545230 

-.548412 

-.224090 

-.224887 

.3 

-.549768 

-.548412 

-.204163 

-.203382 

.4 

-.546987 

-.548412 

-.244610 

-.243848 

.5 

-.551681 

-.548412 

-.396033 

-.395210 

.6 

-.546987 

-.548412 

-.244610 

-.243848 

.7 

-.549768 

-.548412 

-.204163 

-.203382 

.8 

-.545230 

-.548412 

-.224090 

-.224887 

.9 

-.548250 

-.548412 

-.235410 

-.235451 

1.0 

-.548395 

-.548412 

-.137099 

-.137103 
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TABLE 5.- CIRCUMFERENTIAL STRESS RESULTANTS FOR A FREELY 
SUPPORTED CYLINDER WITH A HEATED GENERATOR 


s/L 

T 2 /EhaT for - 

n = 0 
(a) 

n = 1 
(a) 

n = 2 

(a) 

n = 3 
(a) 

n = 4 
(a) 

0 

-0.316899 

-.311717 

-0.316912 

-.311717 

-0.126781 

-.124687 

-0.063403 

-.062344 

-0.037307 

-.036674 

0.1 

-0.024860 

-.024696 

-0.025348 

-.025184 

-0.010746 

-.010682 

-0.005923 

-.005891 

-0.004000 

-.003982 

0.2 

0.020494 

.020413 

0.030190 

.020108 

0.007665 

.007630 

0.003399 

.003381 

0.001511 

.001499 

0.3 

0.005148 

.005102 

0.005147 

.005099 

0.002004 

.001982 

0.000846 

.000833 

0.000203 

.000194 

0.4 

-0.000788 

-.000792 

-0.000752 

-.000758 

-0.000308 

-.000313 

-0.000267 

-.000272 

-0.000422 

-.000426 

0.5 

-0.001020 

-.001015 

-0.001007 

-.001004 

-0.000435 

-.000437 

-0.000349 

-.000352 

-0.000484 

-.000488 

H 

-0.000788 

-.000792 

-0.000752 

-.000758 

-0.000308 

-.000313 

-0.000267 

-.000272 

-0.000422 

-.000426 


0.005148 

.005102 

0.005147 

.005099 

0.002004 

.001982 

0.000846 

.000833 

0.000203 

.000194 

0.8 

0.020494 

.020413 

0.020190 

.020108 

0.007665 

.007630 

0.003399 

.003381 

0.001511 

.001499 

0.9 

-0.024860 

-.024696 

-0.025348 

-.025184 

-0.010746 

-.010682 

-0.005923 

-.005891 

-0.004000 

-.003982 

1.0 

-0.316899 

-.311717 

-0.316912 

-.311717 

-0.126781 

-.124687 

-0.063403 

-.062344 

-0.037307 

-.036674 


a Top line - present solution 


Bottom line - analytical solution (appendix C) 
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TABLE 6.- MERIDIONAL STRESS RESULTANT FOR A 
CLAMPED ORTHOTROPIC CONE 


s/ L 

T l/ E l ha l T max 

Present 

solution 

SADAOS 
(ref. 6) 

SALORS 
(ref. 5) 

0 

67.9152 

68.5935 

71.52 

.i 

34.8894 

34.9249 

34.90 

.2 

23.1281 

23.2835 

23.26 

.3 

17.5112 

17.4625 

17.45 

.4 

13.9475 

13.9699 

13.96 

.5 

11.6352 

11.6416 

11.63 

.6 

9.98591 

9.97855 

9.969 

.7 

8.72492 

8.73126 

8.727 

.8 

7.76298 

7.76115 

7.757 

.9 

6.98650 

6.98525 

6.981 

1.0 

6.49627 

6.49742 

6.680 


TABLE 7.- STRESS RESULTANTS IN A CLAMPED 
TRUNCATED HEMISPHERE 


s/L 

Tj/Eh 

aT 

max 

T 2/ Eh “ T max 

Present 

solution 

SALORS 
(ref. 5) 

Present 

solution 

SALORS 
(ref. 5) 

0 

-0.04843 

-0.05115 

-0.5145 

-0.5155 

.1 

-.05076 

-.04968 

-.1583 

-.1571 

.2 

-.05077 

-.05065 

.06259 

.06220 

.3 

-.05339 

-.05445 

.07233 

.07230 

.4 

-.05813 

-.05895 

.06007 

.06046 

.5 

-.06430 

1 

.06168 

.06225 

.6 

-.07224 


.07493 

.07550 

.7 

-.08469 

-.08630 

.11291 

.1129 

.8 

-.10697 

-.1065 

.12873 

.1272 

.9 

-.12065 

-.1185 

-.22234 

-.2239 

1.0 

-.03444 

-.03821 

-1.0103 

-1.012 
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TABLE 8.- MOMENT RESULTANTS FOR A CLAMPED ANNULAR PLATE 


s/L 

Mi /Eh 3 aT 
L max 

M«/Eh 3 aT 
2/ max 

Present 

solution 

Analytical 
solution 
(appendix C) 

Present 

solution 

Analytical 
solution 
(appendix C) 

0 

1.80057 

1.80055 

2.62350 

2.62350 

.1 

3.11838 

3.16188 

3.39046 

3.53529 

.2 

4.28310 

4.33328 

4.20909 

4.37560 

.3 

5.33251 

5.37015 

5.06365 

5.18880 

.4 

6.29431 

6.31184 

5.94524 

6.00328 

.5 

7.18740 

7.18419 

6.84879 

6.83769 

.6 

8.02591 

8.00600 

7.77156 

7.70467 

.7 

8.82060 

8.79099 

8.71221 

8.61293 

.8 

9.57974 

9.54944 

9.67016 

9.56858 

.9 

10.3098 

10.2892 

10.6454 

10.5760 

1.0 

11.0161 

11.0163 

11.6381 

11.6382 
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TABLE 9.- ELEMENTS OF MATRIX [t r ] 


1 

'* L . ikgiji') 

o 

5c k 

4 

o 


1_ 

€ k L c k R l,k+l\ 

g 

5t k 

4 


f 2 

A 

2 

32R l,k\ 2R l,k/ 


32 

64R l,k 


64 

2 

32R l,k+l\ l 2R l,k+l/ 


32 

64R l,k+l 


64 

15 

1 (-7 ♦ e “ R ‘.A 

o 

7 

€ k 

o 

e k 

15 

1 L ) e k R l,k+l\ 

r» 

_ _7_ 

e k 


E k 

'*k 

I6R l,k\ 2 R l,k/ 


' 16 

32R l,k 


32 

8e k 

16R l,k + l\ 2R l,k+l J 

w 

16 

32R l,k+l 


32 

0 

i L , 

o 

3 

1 

o 

1 

o 

1 ( 3 , e k R l,k+l^ 

ft 

. 3 



1 


4e k R l,k\ 2 R l,ky 


4e k 

~ 8R l,k 


" 8 


4t k R l,k+l\ 2R l,k+l/ 

U 

4£ k 

8R l,k+l 


" 8 

5 

1 4 - 

0 

5 

1 

o 

1 

5 

i / 5 ( e k R i,k+i\ 

g 

5 

1 


1_ 

e 3 

£ k 

2e k R i, k \ 2R i *) 


H 

4e k R l,k 


4e k 

' £ S 

2 R lMl) 


H 

4e k R l,k+l 


4e k 

0 

l L e k R i,k\ 

o 

1 

1 

o 

1 

o 

1 L , £ k R l,k+l\ 

g 

i 

1 

0 

1 


2e k R l,k\ 2Rl ' k / 


2 4 

4£ i R i,k 


4 e 2 

4e k 


2e k R l f k+l\* 2R l,k+l / 


X 

4£2R i,k + i 


*4 

6 

L ( 3 , £ k R Lk\ 

o 

3 

1 

o 

1 

6 

1 L ] € k R l,k+l\ 

g 

3 

1 

o 


' £ I 

^.kl 2Rr J 


'4 

2 4 R i, k 


" 2e k 

e 5 
e k 

e k R l,k+l \ 2R l,k+l / 


~4 

2e k R l,k+l 



0 

1 

2 

0 

0 

£ k 

T 

0 

0 

0 

1 

2 

0 

0 

e k 

8 

0 

0 

0 

3 

0 

0 

i 

0 

0 

0 

3 

0 

0 

1 

0 

0 


2e k 



4 




2e k 



*4 



0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0 





X 







2£ k 



0 

2 

71 

0 

0 

1 

TS 

0 

0 

0 

2 

%3 

0 

0 

1 

e 2 

0 

0 


e k 



e k 




e k 



e k 



0 

0 

1 

0 

0 

e k 

0 

0 

0 

1 

0 

0 

_ik 

0 



2 



8 




2 



8 


0 

0 

3 

0 

0 

1 

0 

0 

0 

3 

0 

0 

1 

0 



X 



"4 




2e k 



4 


0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 






X 







2e k 


0 

0 

2 

4 

0 

0 

1 

4 

0 

0 

0 

2 

'4 

0 

0 

1 

4 
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(a) Coordinate system and displacement component directions. 



(b) Detail of shell wall. 

Figure 1.- Geometry of a shell of revolution. 
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r * 1.0 m 
L * 1.0 m 
fi * 0.25 



E * 10 GPa 
a m 1 x 10“ 5 /K 
h * 10 cm 


V T i- 


sinh (s/L) 

sinh (1/2) (° = S /L < 1/2) 


sinh (1-s/L) 
sinh (1/2) 


(1/2 s s/L 5 1) 


(a) Cylinder with heated midlength - clamped edges (cylinder I). 



h * 1 cm L = 1.0 m a * 1 x 10" 6 /K 

r * 1.0 m E * 10 GPa ji * 0.3 

T q * TV* £ T cos n0 

— - tanh 77 (n * 0) n 

o 

tanh it (n > 0) 

^(1 +n ) 



(b) Cylinder with heated generator - freely supported edges (cylinder n). 



h « 0,0625 cm 
/3m 30° 


L : 
r 


120 cm 
60 cm 


mm 
Ej « 50 GPa 

E 0 * 40 GPa 


T 0 = 1500 - 1000 (s/L) + 500 (s/L) 2 
TV « 500 

fij * 0.25 

^2 ■ 0.20 

£*1 * 6.67 x 10‘ 6 /K 
a 2 “ 6.08 x 10- 6 /K 


(c) Conical frustum — clamped edges. 



h * 1 cm 
L * 1.0 m 
a « 1.0 m 
E * 10 GPa 


a = 1 x 10‘ 6 /K 
H x 0.3 
R x 0.540 m 


t T » T-* 1 + 0.6587 log 1 ± sin (s / L > 1 + C0S (s/L) 

r a ^ o 1 1- sin (s/L) 


j_ 1 - sin (s/L) J 

(d) Truncated hemisphere - clamped edges. 



a = 1.0 m 
b = 2.0 m 
h » 1 cm 
E * 100 GPa 


a x 1 x 10' 3 /K 

L x r 2 - i-i 

fJL * 0.3 

• 2000 (1 + s/L) 2 


(e) Annular plate - clamped edges. 

Figure 2.- Models used for sample calculations. 
All temperatures are in K. 
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Dimensionless stress resultants, Tj/EhaT and T 2 /EhaT 


i 



Figure 3.- Stress resultants in cylinder with prescribed midlength 
temperature. Both edges clamped; cylinder I. 
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Figure 4.- Circumferential stress resultant in cylinder with heated generator. 

Both edges clamped. 



Dimensionless circumferential stress resultant 







Dimensionless moment resultants 


i 

i 

i 



0 .2 .4 .6 .8 1.0 

s/(b-a) 


Figure 8.- Moment resultants in annular plate. Both edges clamped; 

/ \2 
T 0 = T i = 2000 (l +j-J . 
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Dimensionless axial stress resultant, T 1 /EhaT 





